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Abstract 

The two body non-leptonic decays are analyzed in factorization 
approximation, using quark model, ^ = as a free parameter. 

It is shown that the experimental branching ratio for A;, — > AJ/ip 
restricts ^ and this ratio can be understood for a value of which lies 
in the range < ^ < 0.5 suggested by two body B meson decays. The 
branching ratios for Af, — > AcD*{Ds) are predicted to be larger than 
the previous estimates. Finally it is pointed that CKM-Wolfenstein 
parameter + rf' , where r/ is CP phase, can be determined from the 
ratio of widths of A;, — > AD and Af, — > AJ /ip or that of Af, — > pDg 
and Afj — > Ac Dg independent of the parameter ^. 

PACS: 13.30. Eg, 11.30.Hv, 12.39.Jh 
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1 INTRODUCTION 



Two body non-leptonic decays of bottom baryons provide useful informa- 
tion for QCD effects in weak decays and indirect CP asymmetries which 
involve CKM-Wolfenstein parameters p and rj. The standard frame work 
to study non-leptonic decays of bottom baryons is provided by an effec- 
tive Hamiltonian approach, which allows a separation between short- and 
long-distance contributions in these decays. The latter involves the matrix 
elements < MB'\Oi\B > at a typical hadronic scale, where Oi is an operator 
in the effective Hamitonian. These matrix elements cannot be calculated at 
present from first principles. Thus one has to resort to some approximate 
schemes. Such schemes are often complicated by competing mechanisms, 
such as factorization, baryon pole terms and W-exchange terms, each of 
which has uncertainties of its own. The purpose of this paper is to study 
a class of two body bottom baryon non-leptonic decays in the framework of 
factorization scheme, where neglecting final state interactions, hadronic ma- 
trix elements are factorized into a product of two matrix elements of the form 
< B'\J^\B > and < 0| J^|M > for which more information may be available. 

Following the phenomenological sucess of factorization in the heavy to 
heavy non-leptonic B-meson decays this frame work has been extended 
to the domain of heavy to light transitions [Q. The factorization anstaz 
here introduces one free parameter, called ^ = ( being number of 
colors), which is introduced to compensate for the neglect of color octet-octet 
contribution in evaluating the hadronic matrix elements in the heavy to light 
sectors. The range < ^ < 0.5 has been found to be consistent with 
data on a number of measured B meson decays. We apply the factorization 
to decays A;, — ^ AJ/V^, — > A^DsiD*), Ab — > AD and A^ — > pDs. 
In addition, we use quark model to fix current coupling constants which 
appear in the matrix elements < B'\Jn\B >. We show that the measured 
branching ratio for A;, — > AJ/ip can be accounted for in this approach with 
the parameter ^ in the above mentioned range. Our estimates for branching 
ratios for A5 — > AcDs{D*) are larger than their previous estimates |]. 
The decays A;, — > AD and A5 — > pDs can give information on the CKM- 
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Wolfenstein parameter (p^ + ^7^) or iKib/Ktl independent of ^. 
We write the effective HamiltonianlBI: 



H^ff{AB = l) = ^ 



E V^tV*, {C,Ot + COD 

'=U,C 



q=u,c 

where Ci are Wilson coefficients evaluated at the renormalization scale /i; the 
current- current operators Oi^2 are 

0^2 = ic''h)v-A{^'l»)y_^ (2) 

and OJ* are obtained through replacing c by u. Here a and /3 are S'?7(3) color 
indices while (c"6/3)y_A = c°'7^( 1 + 75)6/3 etc. The related Wilson coefficients 
at /i = 2.5GeV^ in next-to-leading logarithmic (NLL) precision are 

Ci = 1.117 

C2 = -0.257 (3) 

These are not very different from those at /i = 5GeV in the leading logarith- 
mic approximation (LLA) 0: Ci(mfe) = 1.11 and C2{mb) = —0.26. 

In the factorization scheme we encounter matrix elements of the form 

<B{p')\J,\B,{p)> = u{p')r,u{p) 

= u{p')i [{gv{s) - ^^(5)75) 7m 

+ {fv{s) + /iA(s)75) (yt^uQu 

+ i{hvis)- fAish5)qf.]u{p) (4) 

where i?b is a baryon, which contains b quark while B is any baryon not 
containing it. Here s = —q^ = —{p — p'Y- In the heavy quark spin symmetry 
limit 1^, the vector and axial vector form factors are related [when B^ belongs 
to the triplet representation of flavor SU(3)] as follows: 

9v{s) = gA{s) = fi (5) 
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fv — hv — hA — —fA — /2 

ruB, 



(6) 



For a decay of the type Bb{p) 
are of the form 

G' 



B{p') + X{px), the matrix elements 
1 I mm' 



T = ^-^ < o| > -(p')r.^p) (^^^ (7) 

In the rest frame of Bh, the decay rate of Bh and its polarization are given 

by 



G' 



/ dsp'{s){p{s)r^{s) + a{s)r'^{s)} 



2 47rm2 

where q — Px — P — p' , s — —q^ and 

r^{s) = {Qis)[g^y + gl)-3mm's{g^y-gl 
+3s {m + m') (^{m — m'Y — gvfv 
- (m - m') ({m + m'f - s) Qa/a] 
+s [Q"{s) (/^ + hi) - Smm's (/^ - h^) 
—2mp'{s)n ■ s — m'^) — 2s) gA9v 

+s {{m - 3m')gvhA - (m + 3m')gAfv) 
+sfvhA (s — m"^ — 5m'^)j | 



Hs) = {Q'is){gl. + gl)-s[im-m'){im + m'y-s)gvfv 

+ (m + m') ^(m — m')^ — 5^^/^ 
1 



(8) 



(9) 



■((m + m'y - s) hi + ((m - m'f - s) fl 
—2mp'{s)n ■ s (m^ — m'^)gvgA 

-s {{m + m')gAhv + (m - m')gvgA) + s^/^-y/A] } (10) 

Here it is understood that form factors are functions of s and p = pv, Pa or 
according as X is 1~, 1"^ or 0~ , while correspondingly a = 0, aA or cr^.g and 

a"' 



= 2^{ 



4m m 



2._./2 



1 



(11) 



Q{s) 
Q'is) 
Q"(s) 



m — m 



) +s (ni^ + m'^) - 2s' 
m — m \ — s [m + m 1 
2 {m^ - m'^)^ - s (m^ + m'^) - s' 



(12) 
(13) 
(14) 



The form factors defined in Eq.(4) are calculated in quark model at s = 
—q"^ = m\ where X is a vector or pseudoscalar particle in the decay B — > 
B'X, thereby taking into account recoil correction. This is in contrast to the 
use of the non-relativistic quark model for the evaluation of form factors at 
zero recoil ^ = [9] . This latter approach also necessitates the extrapolation 
of form factors from maximum q^[—q^ — tm — {rriB — mB'Y] to the desired 
s — —q^ — rn?x. We may point out that since |q| f» l.lbGeV in — > 
A -\- J/ip, for example, the no recoil approximation does not seem to be 
justified; in fact |q| ^ in A, making the s quark in A relativistic. In 
our approach no recoil approximation, nor any extrapolation of form factors 
at the physical point arc needed. Our quark model results do satisfy the 
constraints imposed by heavy quark spin symmetry. 

The plan of the paper is as follows: Sec. 11 summarizes the calculation 
of the baryonic form factors within the frame work of quark model at the 
desired value oi s — —q'^, rather than at the zero recoil point, relegating the 
details in the appendix. In section III we apply the results to some specific 
nonleptonic decay modes of A^. Section IV summarizes our conclusions. 



2 BARYONIC FORM FACTORS IN QUARK 
MODEL 

In order to calculate the form factors we first reduce the matrix elements in 
Eq.(4) from four component Dirac spinors to Pauli spinors without making 
any approximation and do the same for the quark level current 

= + 75)& (15) 



5 



We treat the b quark in Bb extremely non relativistically {pb/rrih ~ 0) and 



set Pb — Pg = q = — p', Eq = y'q^ + m^, E^ = rrib = m^. Then as shown in 
the appendix: 

hA = -/a, fv = hv (16) 

gv{s) = CvIa{E',E',) 

Ms) = -^vlKE',E',) 
m 

gA{s) = UIa{E',E',) 

Us) = --UIb{E',E',) (17) 
m 



where 



a{E\E',) 



1 fW{E' + m'){l-^) + {E', + m',){l + ^) 



2V^3 ^{E' + m'){E', + m',) 



HE' E') ^ - [E' {E' + m')-{E', + m',) 
' 2V£^^ y'(E' + mO(i?^ + m^) 

and E' — p'g, E'^ — Eg — ^J-p''^ + , and mg = m^. Note the exphcit 
appearance of 1/m corrections in the above formulae. Here and are 
respectively the spin-unitary spin part of the matrix elements of the current 
operator (13); for example, for Bi, belonging to the triplet representation of 
<S't/(3), = iA and / is overlap integral 

I = N^Ni I rf (pi2, k - ""'^r' p') ^»(Pi2> ^)d'pi2d'k (19) 

The recoil correction is represented by momentum mismatch "^^tj"^ p^, which 
arises since the rest frame of Bj, is not that of Bg baryon. Here m' = mi + 

m2 + mg where mi and m2 are masses of the spectator quarks and m'^ is 
that of q quark resulting from the decay of b. Note that the form factors in 
Eq.(15) are determined at the desired value of s = — g^. 

As already noted for B^ belonging to the triplet representation = 
and then the relations (14) and (15) are consistent with those given in Eqs.(5) 
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Table 1: Quark model predictions for baryonic form factors for A;, transitions. 
P = 0.51 GeV, P' = 0.44 GeV for p and A and = 0.48 GeV for Ac. h^gv^ 



9a, /2//1 = 


ifv/gy) 


m = {hy/gv)m = 


{hA/gA, 


m = — 




m. Note that 


only the last column 


depends 


on the overlap integral /. 






Transition 


\p'\ 




fi{s)/I 


/2//1 


^1 




/ 




pDs 


2.376 


I/V2 


0.720 


0.123 


^ 1 




0.119 


0.086 


AD 


2.374 


-1/Vs 


-0.558 


0.129 


1 


1 


0.215 


-0.120 




1.756 


-i/Vs 


-0.604 


0.158 


0.943 


0.826 


0.426 


- 0.257 




1.766 


1 


1.052 


0.131 


0.978 


0.983 


0.791 


0.829 




1.850 


1 


1.048 


0.137 


0.949 


0.908 


0.810 


0.852 



and (6) obtained in the heavy quark spin symmetry limit. To proceed further 
we use harmonic oscillator or Gaussain wave functions in Eq.(20)to obtain 



2(3(3' 




■4 m'2 2(/32 + /5'2) ^ 



+ 

We take P or P' as [10] 

P' = (21) 

where /iq — J^^^^^ is the reduced mass of the bound system, M^v being 
the nucleon mass and Mh that of B, D, K* or p meson for A?,, Ac, A and p 
respectively, k is the spring constant and its value is taken to be (440 MeV)^ 
[11]. 

We summarize in Table 1, the form factors gvi^) — gA{s) — /i, fv{s) — 
hv{s) = /lyi(s) = — /a(-s) = f2/m for the transitions A^ — > pDs, A}, — > AD, 
Aft — > AJ/^, Afe — > AcD;{Ds), for s = mjy^, m^j/^, and (m^J, 
"^3 — "^u, rris and rric respectively. For the numerical work we have taken the 
relevant masses ( in GeV) as m = uia^ — 5.641, uia — 1.1157, niA^ — 2.285, 
rup = 0.938, mj/^ = 3.097, mo* = 2.112, mo = 1-864, mo, = 1-968, 
nis — 0.510, nic — 1.6 and m„ = 0.340. 
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3 APPLICATIONS 



We consider those decays of A;, for which baryon poles either do not con- 
tribute or their contribution is highly suppressed due to Okubo-Zweig-Iizuka 

(OZI) rule and that it scales as inverse of m\^. 

For decays of type Ab(p) — > Bq{p')V{q), where is a vector meson, 

PAis) = - aAis) (22) 

Pv{s)^F^5{s-ml) (23) 

where 

<0|j;|y >=Fve^ (24) 
Then Eqs.(9) and (10), on using the relations (5), give the decay rate 

^ ^ ^^^i|^QW)[2/r (25) 

while the asymmetry 

_ -2m\p'\ [(m^ - m'^) - 2m^] J^^iK) 

2Q{m'y) TX{ml) ^^^^ 

where 

T,[my)-^ 3- A+^g(<y:^/ ^^^^ 

^ ^ \ m m'^ — m'^ — 2my fi m? m? — m'^ — 2my J 

The prediction for a is independent of the value of the overlap integral 
and provide a test of the predictions (14) and (15) with ^y = through the 
presence of /2//i- The corrections due to form factors which scales as 1/m 
are dumped into JF functions 

If the vector meson V is replaced by a pseudoscalar meson P, then 

py{s) = = pa{s) (29) 
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cr^(s) = Fl5{s - ml) (30) 

where 

< 0\f^\p >= Fpq, (31) 
Then Eqs.(9) and (11), on using the relations (5), give 

= ^^P^Q'(^pWiiml)]J'[{ml) (32) 
^ -2m\p'\[{m'-m'^)]J^[iml) 



where 
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^1 {^p) = |1 - — 'qI^) + ^ 2Q'{ml) Jl] 

(34) 



2m' mp f2 mp f. 



We are now ready to consider the specific decays. We first consider — ^ 
KJ/if)^ where the first part of the Hamiltonian (1) with q = c and Fierz 
rearrangement give 

G' = GFV,tK:{C2 + ^Ci) (36) 
j; = c7^(l + 75)c (37) 
The constant Fj^^ is determined from T{J/ip — > e+e~) = (5.26±0.37) KeV 



= 1.637 X 10-1 GeV'^ (38) 

Using Gf = 1.16639 x 10~^GeV~^ and = 0.0393 ± 0.0028, \V^s\ = 
1.01 ± 0.18, we obtain from Eqs.(23) and (24) 

r = 8.21xlO-'\C2 + ^G,fmm]^^)J^rim]/^) (39) 
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This gives the branching ratio 



B (A, AJ/^) = 1.47 X 10-^ {C, + iCf fl{m],AT'({m],A (41) 



where we have used Q Ta, = 0.847 x lO^^s"^ = 5.59 x IQ-^^GeV. Using 
Table 1 we finally obtain 

B (Afc — > AJ/iP) = 9.14 X 10"^ (C2 + ^Cif (42) 

a {Ab — > AJ/tp) = -0.18 (43) 

In Fig.l, we show the branching ratio B{Ah — > AJ/ip) as a function of 
^. This decay mode is sensitive to ^ and comparison with the experimental 
value [13] (3.7 ± 2.4) x 10"^ shows that ^ is restricted to < ^ < 0.125 or 
0.35 < ^ < 0.45, which lie within the range < ^ < 0.5 suggested by the 
combined analysis of the present CLEO data on B — > hih2 decay 0. We 
may remark that /2//1 correction to the decay rate is about 6% while that 
to the asymmetry parameter a is about 14%. 

Other decays of interest for which the first part of Hamiltonian (|I]) with 
g = c is responsible are A;, — > ^t^s A;, — > A^D*~. For these decays 

G' = GFV^bK: (Ci + eCa) (44) 

and 

j; = s7;.(l + 75)c (45) 

Then Eqs.(23) and (24) and (29) and (30) [on using the relations (|)] give 
respectively 

r(A,^A+Z}:-) = 2.12 X 10-'' (G, + ^G^)' fK^lt)^riml,) 

(46) 

{Ab^AtDl-) ^ -0.42^^ (47) 
r(A,^A+D;) = 1.50xW-''{G,+^G2ffKmlJJ'[{mlJ 

(48) 



a 
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Table 2: Predictions for the branching ratios (BR) in % for — > ^t-^s 



and Afe — > in the large limit = 0). 



Decay processes 


Present BR calculation 

(e = o) 


BR 
Ref.i 


BR 
Ref.i 


Afc AtD*- 
A, AtD; 


3.26 
2.23 


1 70+0.20 

-L- ' -J-O-SO 

Qn+0-30 

^•'JU_o.40 


1.77 
1.156 



a (A, ^ AM = -0.98 Jp) (49) 

Here we have used F^^ = F^* = 2?>2 MeV\\.2] (in the normalization F^^ = 
131MeV"). Using Table 1, the above equations give 



B (Afc — > A,Dl) = 2.61 (Ci + ^C2f x 10"^ (50) 

a (Afc — > A^Dl) = -0.40 (51) 

BiAb^A.D,) = 1.79 (Ci + ^^2)' X 10-2 (52) 

a (Ab — > A^D,) = -0.98 (53) 

The above branching ratios are not sensitive to ^: 2.55 x 10^^ < B{D*) < 
3.26 X 10-2 and 1.75 x 10^^ < ^p^) < 2.23 x 10"^ for 0.5 > ^ > 0. The /2//1 
corrections are negligible when the meson in the final state is while for 
1^ they are about 5% for the decay rate and for the asymmetry parameter a. 
Previously the above decays have been analyzed in the heavy quark effective 
theory (HQET) with the factorization approximation in the large limit 
either by parameterising the Isgure-Wise form factor Gi{v ■ v') [c.f. Eq.(|^) 
with /i = Gi + {rn/^^/m\JG2, f2 = — G'2/mAj, where since Ac, Af, form a 
multiplet, the absence of second class currents implies G2 = 0] or by 
evaluating it in the large Nc limit |Q. In contrast we have used quark model 
to fix the baryonic form factors and Eqs.(15) and (16). The comparison of 
our predicted results with the previous results mentioned above is presented 
in Table 2. 

Finally we consider the decays Af, — > AD° and A;, — > pDs\ the in- 
terest here is that the ratio of their decay widths with A;, — > AJ/ip and 
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Afe — > Ac-Ds respectively can fix the CKM-Wolfenstein parameter (p^ + rj'^) 
or |14fe/Kf)|, independent of ^, wfiere rj indirectly determines CP- violation. 
For these decays the second part of the Hamiltonian([T|) with q = c [and Fierz 
rearrangement for the former] give 



r(A„ 



AD' 



Gf 



V2 



fr {ml)f:F[{ml)Q'{ml 



(54) 



r(Ab 



2|p'| ^ 



47rml 



(55) 



Using Table 1, Fd = 200 MeV and taking into consideration differences in 
phase space factors p', Q and Q' we obtain 



r A 



AD° 



r (A, A j/^) 

r(A, — .pD,) 
r(A6 ^ A,D,) 



5.88 X 10^ 



ub 



2 X 10~2 



c6 



2.8 X 10"^(p^ + r7^) (56) 
= 9.7x 10-^(p2 + r/2) (57) 



4 CONCLUSIONS 

We have analyzed some two body non-leptonic A;, decays in the factorization 
approximation, treating ^ = l/Nc (which is supposed to compensate for the 
neglect of color octet-octet contribution in evaluating the hadronic matrix 
elements) as a free parameter. In addition we have used quark model to fix 
baryonic form factors at the desired value of s = — without making no 
recoil approximation. The form factors obtained are consistent with the pre- 
dictions of heavy quark symmetry and explicitly display 1/mf, corrections. 
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The experimental branching ratio for — > AJ/ip restricts ^ and can be 
understood for either < ^ < 0.125 or 0.3 < ^ < 0.45. Our predictions 
for the branching ratios — > AcDg{D*) are larger than the previous esti- 
mates. Future experimental data from colliders are expected to verify and 
distinguish the various results. Finally the parameter | Kb/Kb | or (p^ + 77^) 
can be determined independently of the parameter ^ from the ratio of de- 
cay widths of Af, — > AD and A;, — > A J/ i]j or that of Af, — > p Dg and 
A5 — > AcDg, although the branching ratios expected for these decays may 
be hard to measure. 

Wc want to emphasize that our derivation of Eqs. (16) and (17) does not 
depend on the details of quark model. The basic assumption is that in the 
heavy quark limit, the velocity of heavy quark can be neglected. The details 
of the quark model enter in the derivation of the overlap integral /. 

It may be noted from the structure of Eqs. (9) and (10), that the contri- 
bution of the form factors fv, hy, Ha and are proportional to ^ [same 
is true for the term containing {gy — g\)]. Hence when ^ ^ 1, their con- 
tribution can be neglected and in this case asymmetry parameter a is given 

by 



a ~ 



9v + 9a 
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A APPENDIX 

We outhne the derivation of relations (16) and (17). We first reduce the ma- 
trix elements in Eq.(4) from four component Dirac-spinors to Pauli spinors. 
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Thus in the rest frame of Bi, 



IE' + m' 



2E' 



< B(p')\J,\B,(p) > 



gv(s) - qohv(s) - ^,^^J v{s) 



+ 



hA{s) + -^7^ {9Ais) - qofAis))\ a • q| {A - 1) 
< B{p')\J\B,{p) > 



IE' + m' 



2E' 
1 + 



-9a{s) +[qo + — 



E' + m' 



hA(s) 



p;^ ; ) fv{s) + , gv{s) 

E' + m' J E' + m' 

+ F/ ! / (9v{s) + qofv{s)) 

Hi -t mi 



E' + m' 



y[/iA(s) + /A(s)]q(7-q 



icr X q 

q 

{A -2) 



where E'{s) = p'o{s), q = — p', qo = \/\^.^ + s- It may be noted that no 
approximation has been made so far. On the other hand, the Pauh reduction 
of the quark level current 



is given by [ with pQ = p^, pg = p'^]: 



(A -3) 



Jo 



2[E-iE',{E-, + ms){E', + m'^)]'/' 
X {(E'^ + m'^)(Es + ma) + p'g -ps + ia- (p'g x pg) 
{E'2 + m'^)a-ps-{E3 + m3)a-p's} (^-4) 



J = 



2[E,E^{E3 + m3){E', + m',)]'/'' 

X {[-(^3 + "^3)(^3 + "^3) + P3 • Ps] (7 + lip's X Ps) 
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-{(^ ■ PaJPa - ■ P3P3 

+{E'^ + m'^) (p3 -ia X pg) 

+ {E3 + m3){p's + iaxp's)} (A - 5) 

We now treat the quark Q extremely non-relativistically and thus put Ipal — 
0. Then 

jo- {{E', + m',) - a ■ p'^} {A-6) 

pE'^{E'^ + m'^) 

J = ^— ===={-(^3 + "^3)^ + P3 + ^^-P3} (^-7) 

where 

E', = ^Jp'l + m'i = V(P3 - q)2 + mf ^q^ + m'i 

Suppose that the initial baryon B contains a heavy quark Q {b in our case) 
and two light quarks qi and q2 which behave as spectators. The final baryon 
B' is composed of the quark q [s, c, or u quark] and the same spectators as in 
B. For the initial baryon composed of quarks Q{= 53), qi, q2, we introduce 
relative coordinates and momenta as 

P12 Pi P2 17111712 

ri2 = ri - r2, = , mi2 = 



mi2 mi 777.2 + 1712 

_ miFi + m2r2 _ 

•tli2 — , ri2 3 — ri2 — r3, 

mi2 

k P12 P3 1713(1111 + 1112) 

P12 = Pi + P2, - = ^ , = Z 

1 1711+1712 . . 

m = mi + 1712 + tns, k = — P12 ; Ps - 8) 

m m 

For the initial baryons, its rest frame is its center of mass frame so that 
Pi + P2 + Pa = which implies P12 = — P3 = k and then 

mi 

Pi = P12 H ■ k 

mi + 777.2 

P2 = -P12 + — k (A - 9) 

7711 + ^2 
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Denoting the relative momenta of quarks in the baryon B' by primes and 
noting that p'^ = pi, p'^ = P2 so that = P12, P'12 = P12, giving p'3 = 
— P12 + p' = — k — q and 

k' = k-^^^i±^p' (A-10) 

m' 

Calhng t/^s the spatial wave function in momentum space and noting that 
when p'3 in Eqs.(A-6) and (A-7) is replaced by — k — q, the linear terms 
in k do not contribute in the spatial integral and as such the right sides of 
Eqs.(A-6) and (A-7) are independent of integration variables k, pi2 and k'. 
The comparison of hadronic matrix elements in Eqs.(A-l) and (A-2) with 
those of Eqs.(A- 6) and (A-7) give the relations (16) and (17). The use of 
delta function 5(pi-p\), 5(p2-p'2), <^(Pi+P2 + P3) and (p'i+p'2+P'3-p') 
reduce the spatial integral to the form given in Eq.(20). 
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Figure Captions 

Figure 1. Branching ratio for A5 — > AJ/ip as a function of ^. The dotted 
and sohd hues show the CDF measurement. 
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